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Abstract. A new proof of Friedrich's theorem on the existence and stabihty of asymptotically 
de Sitter spaces in 3+1 dimensions is given, which extends to all even dimensions. In addition we 
characterize the possible limits of spaces which are globally asymptotically de Sitter, to the past 
and future. 
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1. Introduction. 

a : 

^ . Consider globally hyperbolic vacuum solutions (M"+^, 17) to the Einstein equations with cosmo- 

logical constant A > 0, so that 

in ; 

(N : (1-1) Ricg -^g + Ag = 0. 

> : 

CN I The simplest solution is (pure) de Sitter space on M""^^ = M x S", with metric 

! 

O ■ (1.2) gas = -dt^ + cosl?{t)gsr.,iy 
00 ■ 

I More generally, let {N^,g]\f) be any compact Riemannian manifold with metric g^ satisfying the 

■ Einstein equation Ricg^ = (n — l)gN- Then the (generalized) de Sitter metric 

"o : (1.3) g^s = -dt^ + cosh\t)gN, 

I < on R X A'" is also a solution of (1.1), with A = n(n — l)/2. 

{3Jq| Let dS~^ be the space of all globally hyperbolic spacetimes {M^^^,g) satisfying (1.1), with a 

' spatially compact Cauchy surface, which are asymptotically de Sitter (dS) to the future, i.e. future 

• ^ . conformally compact in the sense of Penrose; the terminology asymptotically simple is also used in 

r> I this context. Thus there is a smooth function Q such that the conformally compactified metric 

^- : (1.4) -g = n-'g, 

extends to the compactified spacetime M = M UX"^, where X+ is a compact n-manifold without 
boundary. The function is smooth on M, with O > 0, = O~^(0) and dU 7^ on 
The boundary metric 7 = g\x+ depends on the choice of fi; however the conformal class [7] of 
7 is independent of and is called future conformal infinity. Such spacetimes are geodesically 
complete to the future of an initial compact Cauchy surface S diffeomorphic to X"*". There are no 
restrictions on the class [7] or the topology of X"*" , and so such spacetimes are sometimes also called 
asymptotically locally de Sitter. Changing the time orientation gives the same notion for dS~ , with 
future conformal infinity X+ replaced by past conformal infinity X~. The smoothness of g in (1.4) 
up to M may be measured in Holder spaces C™'", but we will mostly use Sobolev spaces which 
are more natural in this context. 

In addition, let dS^ be the space of such globally hyperbolic spacetimes which are in both dS'^ 
and dS~; thus such spacetimes are (completely) global, in the sense that they are geodesically 
complete and asymptotically simple both to the past and to the future. 
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Mathematically, the most significant result on the structure of such spacctimcs is Friedrich's 
theorem [8], [9] that in 3 + 1 dimensions, the Cauchy problem with data on X+, (or T~) is well- 
posed, cf. also [10] for recent discussions. Thus, for arbitrary Cauchy data on X+, there is a 
unique spacetime (M^,p) which realizes this data at future infinity. Moreover, small but arbitrary 
variations of the Cauchy data give rise to small perturbations of the solution. It follows in particular 
that the space dS^ of global solutions is open; thus spaces in dS^, in particular pure de Sitter 
space (M^, gds)^ are stable under small perturbations of the Cauchy data at X+, (or T~). The same 
statement holds for perturbations of the data on a compact Cauchy surface S for (M^,^). 

The purpose of this paper is to extend Friedrich's theorem to arbitrary even dimensions. Let X+ 
be any closed n-manifold, n odd, and let 7 be any smooth Rieniannian metric on X"*", s > ^+1. 
Next, let r be any symmetric bilinear form on X+ satisfying the constraints 

(1.5) trjT = 0, 6jT = 0, 

i.e. r is transverse-traceless with respect to 7. Define 71 ~ 72 and ti ^ T2 if these data are 
conformally related, i.e. there exists A : X+ such that 72 = A^7i and T2 = /(A)ri, where / is 

chosen so that (1.5) holds for T2, cf. [5] for the exact transformation formula. Let ([7], [r]) be the 
equivalence class of (7,t). Then Cauchy data for the Einstein equations (1.1) with A > consist 
of triples (X"*", [7], [r]). The form r corresponds to the order n behavior of the metric; roughly for 
g as in (1.3), r = (9q)"'^|j+; see §2 for further details. 

Theorem 1.1. The Cauchy problem for the Einstein equations with Cauchy data (X"*", [7], [r]) at 
future conformxil infinity is well-posed in H^^"^ x H'^ , for any s > ^ + 2. 

Thus, given any Cauchy data ([7],[t]) G H^~^"'{I~^) x H^{Z~^) satisfying (1.5), up to isometry 
there is a unique Einstein metric {M'^'^^^g) G dS^ whose conformal compactification as in (1.4) 
induces the given data ([7], [r]) on X+. 

This result has the following simple consequence: 

Theorem 1.2. The space dS^ is open with respect to the H^^^ x topology on , s > § + 2. 
Thus, given any dS solution (M"^^, (70) G dS^, any H^^"' x H"^ small perturbation of the Cauchy 
data ([7],[t]) on X+ (or X~) gives rise to complete solution {M^~^^,g) G dS^ globally close to 
{M'^'^^^go). In particular, the even- dimensional pure de Sitter spaces g^g in (1-3) are globally 
stable. 

Here, globally close is taken with respect to a natural topology on the conformal compact- 
ification M = M U X+ UX^, sec the proof for details. The complete solution (M"+^,(/) induces 
jjs+n ^ jjs Cauchy data at both past and future conformal infinity X^, X+. Of course the size of 
the allowable perturbations in Theorem 1.2 depends on (M"+^,5o)- 

We describe briefly the main ideas in the proof of Theorem 1.1; full details are given in §2. 
The Einstein equations (1.1) induce a 2"^*^ order system of equations for a compactificd metric g 
in (1.4). However, this system is degenerate at X"*" = {0 = 0} and this degeneracy causes severe 
problems in trying to prove the well-posedness of the system. In 3 -|- 1 dimensions, Friedrich [8] 
has developed a larger and more complicated system of evolution equations, the conformal Einstein 
equations, for the (unphysical) metric g together with other variables. This expanded system is 
non-degenerate and shown to be symmetric hyperbolic; then standard results on such systems lead 
to the well-posedness of the conformal field equations. However, it seems very unlikely that this 
method could succeed in higher dimensions, cf. [10], due at least in part to the special form of the 
Bianchi equations in 3 -f 1 dimensions. 

The approach taken here is to replace the Einstein equation by a more complicated but con- 
formally invariant higher order equation for the metric alone, whose solutions include the vacuum 
Einstein metrics (with A term). In 3 -|- 1 dimensions, this system is the system of 4*^ order Bach 
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equations, cf. (2.12) below. The Bach equations have been used in a number of contexts in con- 
nection with issues related to conformal infinity, cf. [14], [15], [16], for example. 

In higher even dimensions, in place of the Bach tensor, we use the ambient obstruction tensor 
H of Fefferman-Graham [6], which agrees with the Bach tensor in 3 + 1 dimensions; this tensor is 
also characterized as the stress-energy tensor of the conformal anomaly, cf. [5] . The tensor 'H is a 
symmetric bilinear form, depending on a given metric g on M""*"^ and its derivatives up to order 
n + 1. The equation 

(1.6) H = Q 

is conformally invariant, and includes all Einstein metrics (of arbitrary signature and A-term). It is 
a system of (n + l)'^*-ordcr equations in the metric, whose leading order term in suitable coordinates 
is of the form □~2~, where □ is the wave operator of the metric g. Conformal invariancc implies 
that the system (1.6) is non-degenerate at X"*" = {Q = 0}. Theorem 1.1 is then proved by showing 
that natural gauge choices for the diffeomorphism and conformal invariance of (1.6) lead again to 
a symmetrizable system of evolution equations. 

In the context of Theorem 1.2, it is of interest to understand the closure dS"^ of the space dS^, 
i.e. the structure of spacetimes which are limits of spacetimes in dS^ but not themselves in dS^. 
A first step in this direction was taken in [2] in 3 -|- 1 dimensions, and Theorem 1.1 allows one to 
extend this to any even dimension. Let dS^ be the closure of dS^ with respect to the x 
topology on the Cauchy data on either or T~, i.e. the union of the closures with respect to data 
on X~ and . Let ddS^ = dS^ \ dS^ be the resulting boundary consisting of limits of spaces in 
dS^ which are not in dS^. 

Theorem 1.3. For (n + 1) even, a space in the boundary ddS^ ofdS^, is described by one of the 
following three configurations: 

I. A pair of solutions {M,g~^) G dS~^ and {M,g~) € dS~ , each geodesically complete and globally 
hyperbolic. One has Z~ = for {M,g~^) and = for {M,g~). Both solutions {M,g~^) and 
{M,g~) are "infinitely far apart". 

II. A single geodesically complete and globally hyperbolic solution {M,g) G dS'^ , either with a 
partial compactification at I~ , or I~ = 0. 

///. A single geodesically complete and globally hyperbolic solution (M, y) G dS~ , either with a 
partial compactification atl'^, orX+ = 0. 

Cases II and III have been distinguished here, but these behaviors become identical under a 
switch of time orientation. 

One of the main points here is that singularities do not form on spaces within dS^. One does 
expect singularities to form "past" the boundary ddS"^. The most natural limits are those of type 
I; this behavior occurs very clearly and explicitly in the family of dS Taub-NUT metrics on R x S^, 
cf. [2] for further discussion. It would be very interesting to know more about the structure of 
dS^; for instance, is it compact and connected? 

Theorems 1.1-1.3 are proved in §2, and we close the paper with some remarks on extending these 
results to vacuum equations with A < and to the Einstein equations coupled to matter fields. 

I would like to thank the referee and Piotr Chrusciel for very useful comments on the paper. 

2. Proofs of the Results. 

Throughout the paper, we consider globally hyperbolic vacuum spacetimes (M, g) with A > in 
(n + 1) dimensions. By rescaling if necessary, it is assumed that A is normalized toA = n(n — l)/2, 
so that the Einstein equations read 



(2.1) 



Ricg = ng. 
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The simplest solution of (2.1) is (pure) deSitter space on M = R x 5", with metric (1.2), or 
its generalization in (1.3). These de Sitter metrics gds are geodesically complete and globally 
conformally compact, i.e. in dS^. In fact, defining s € (—§,§) by cosh(t) = 3^^^ and letting 

(2.2) g = cos2(s)5, 
one has 

9ds = -ds^ + gs^{i), 

which is the metric on the Einstein static spacetime in the region s G [— f , f ]• The metric g^s is 
real analytic on the closure M = M U J+ U X~ and the loci T"*" = {s = |} = {t = 00}, X~ = {s = 
— ^} = {t = —00} represent future and past conformal infinity. The induced metric on is of 
course the unit round metric gs'^(i) on S"^. The same discussion holds for {N,gi\f) as in (1.3) in 
place of gsn{i). 

Consider Einstein metrics (M"+^, (7) in dS^, so that there is a compactification 

(2.3) g = p^g 

as in (1.4) to future conformal infinity with = {p = 0}; all of the analysis below works 
equally well for spaces in dS~ . 

A compactification g = p^g as in (2.3) is called geodesic if p{x) = distg{x,I~^). These arc often 
the simplest compactifications to work with for computational purposes. Each choice of boundary 
metric 7 G [7] on X"*" determines a unique geodesic defining function p, (and vice versa) . The Gauss 
Lemma gives the splitting 

(2.4) g = -dp'' + gp, g = p'^-dp"" + gp), 

where gp is a curve of metrics on The asymptotic behavior of g at X+ is thus determined by 
the behavior of gp as p ^ 0. For example, the geodesic compactification of the de Sitter metric 
(1.2) with respect to the unit round metric at 1^ is 

gds = -dp'' + {l + {^ffgsn(i), 

for p G [0, 00). 

Now consider a Taylor scries type expansion for the curve gp on This was analysed in case of 
asymptotically hyperbolic or AdS metrics with A < by Fefferman-Graham [6], and for dS metrics 
by Starobinsky [19] when n = 3. This idea of course has further antecedents in the Bondi-Sachs 
expansion and peeling properties of the Weyl tensor when A = 0. In any case, the FG expansion 
holds equally well for metrics in dS^ (or dS^) in place of asymptotically AdS metrics; in fact the 
two expansions are very closely related, cf. [2], [18] and further references therein. 

The exact form of the expansion depends on whether n is odd or even. If n is odd, then 

(2.5) gp ~ 5(0) + p^5(2) + -.■ + P"~^5(n-l) + P"5(n) + P"+^5(n+l) + -., 

with 5(0) = 7. 

This expansion is even in powers of p up to order n — 1. The coefficients (7(2fc), < k < n/2 arc 
locally determined by the boundary metric 7 = 5(0); they are explicitly computable expressions in 
the curvature of 7 and its covariant derivatives. For example for n > 3, 

(2.6) ,p, = _i_(ffi,,-_A_,), 

cf. also [5], [2] for formulas for ^(fc) for A; > 2. 
The term 5(„) is transverse-traceless, i.e. 

(2.7) trj9{n) = 0, <^75(n) = 0, 

4 



but is otherwise undetermined by 7 and the Einstein equations (2.1); thus, at least formally, it is 
freely specifiable. For k > n, terms ^(fc) occur for k both even and odd; the term g(^i^^ depends on two 
boundary derivatives of g(k-2) ■ The main point is that all coefficients g^^^ are locally computable 
expressions in ^(q) and . 

Mathematically, the expansion (2.5) is formal, obtained by compactifiying the Einstein equations 
and taking iterated Lie derivatives of ^ at p = 0. If the geodesic compactification g is in C""'"(M), 
then the expansion holds up to order m + a, in the sense that 

(2.8) gp = 3(0) + p'5(2) + •••• + P™3M + 0(p™+°). 
Suppose instead n is even. Then the expansion reads 

(2.9) gp - 5(0) + P^9{2) + ■■■■ + p""^5(n-2) + p"5(n) + P"(logp)7^ + ... 

Again the terms 5(2^) to order n — 2 are explicitly computable from the boundary metric 7, as 
is the coefficient V. of the p^i^ogp) term. The term g(y^ satisfies 

tr^a{n) = a, 6jg(^n) = b, 

where a and b are explicitly determined by the boundary metric 7 and its derivatives, but g'(„) is 
otherwise undetermined by 7 and the Einstein equations; as before, it is formally freely specifiable. 
The series (2.9) is even in powers of p, (at all orders) and terms of the form p^^(\ogpy appear at 
order > n. Again the coefficients g^j^-^ and Tif^k) depend on two derivatives of g(k-2) a-^d 7^(fc_2). 

Although the expressions (2.5) and (2.9) are only formal in general, Fefferman-Graham [6] showed 
that if the undetermined terms (5(0)) 5(n)) analytic on the boundary T"^, then the expansion (2.5) 
converges, (for n odd), cf. also [2]. Thus gp is analytic in p for p small and one has a dS Einstein 
metric in this region given by (2.4). A similar result has recently been proved by Kichenassamy 
[13], (cf. also [17]), for n even; in this case the polyhomogeneous expansion (2.9) converges to gp 
for p small. 

The term TC, which appears only when n is even, has a number of important interpretations. 
First Fefferman-Graham [6] observed that this tensor, locally computable in terms of the boundary 
metric 7, is a conformal invariant of 7 and is (by definition) an obstruction to the existence of a 
formal power series expansion of the compactificd Einstein metric; in fact it is the only obstruction. 
The tensor 7i is also important in the (A)dS/CFT correspondence, in that (up a constant) it equals 
the stress-energy tensor (i.e. the metric variation) of the conformal anomaly of the corresponding 
CFT, cf. [5]. It also arises as the stress-energy or metric variation of the Q-curvature of the 
boundary metric 7, cf. [7]. 

The tensor TC is transverse-traceless 

(2.10) 6n = trH = 0, 

and a conformal invariant of weight 2 — n, i.e. iig = X^g, then H = X^-'^H. Further, if g is 
conformal to an Einstein metric, with any value of A, then 

(2.11) H = 0. 

In addition, as observed in [6] , these properties hold for metrics of any signature and so the equation 
(2.11) can be viewed as a conformally invariant version of the Einstein equations with an arbitrary A 
term and arbitrary signature. We are not aware of any analogue of such a tensor in odd dimensions. 

We will use the tensor Ti to study de Sitter type solutions of the Einstein equations (2.1). Al- 
though the derivation of the obstruction tensor H. arises from the structure at infinity of conformally 
compactified odd-dimensional Einstein metrics, once it is given, one can use it to study the Einstein 
equations themselves in even dimensions. Thus, replacing ra by ra -|- 1, a vacuum solution of the 
Einstein equations (M"+^,5f) with A > 0, (or any A), in even dimensions is a solution of (2.11). For 
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G dS'^, the equation (2.11), being conformahy invariant, also holds for the compactified 
Einstein metric g in (2.2); moreover it has the important advantage of being a non-degenerate 
system of equations in g. As is well-known [8] , the translation of the Einstein equations for (M, g) 
to the compactified setting g leads to a degenerate system of equations for g. 

When n = 3, so dimM = 4, up to a constant factor H is the Bach tensor B, given by 

(2.12) B = D*D{Ric - ^g) + D'^{tr{Ric - ^g)) + 7^, 

where 7?. is a term quadratic in the full curvature of g. (The specific form of TZ will not be of 
concern here). In general, for n > 3 odd, one has, again up to a constant factor, 

(2.13) H = {D*D)'^-^[D*D{P) + D'^{trP)] + L(D"7), 
where 

(2.14) P = P(7) = Ric, - ^g, 

cf. [7] for example. This is a system of PDE's in the metric g, of order n + 1; L(D'"'g) denotes lower 
order terms involving the metric up to order n. 

The Bach equation B = was originally developed by Bach as a conformally invariant version 
of the Einstein equations (with A = 0), and has been extensively studied in this context, cf. [14], 
[15], [16] for some recent work and references therein. It was also used in [3] to study regularity 
properties of conformally compact Riemannian Einstein metrics. 

While Einstein metrics, (of any signature and A), are solutions of (2.11), of course not all solutions 
of (2.11) are Einstein. In addition, for Lorcntzian metrics, 7i is not a hyperbolic system of PDE's 
in any of the usual senses; the equation (2.11) is invariant under diffeomorphisms and conformal 
changes of the metric, and so requires at least a choice of diffeomorphism and conformal gauge to 
obtain a hyperbolic system. 

To describe these gauge choices, suppose {M,g) € dS~^, so that g is an Einstein metric, satisfying 
(2.1), and so (2.11), which is asymptotically dS to the future. Assume that {M,g) has a geodesic 
compactification which is at least C"; then 

(2.15) g = p-^g = -dp'^+gp, 

and gp has the expansion (2.8), with m = n, a = 0. One has = {p = 0} and we set 7 = g(Qy 
By the solution to the Yamabe problem, one may assume without loss of generality that the 
representative 7 G [7] has constant scalar curvature, i.e. 

(2.16) Rry = const, 

on T+. However, closer study shows that the operator P in (2.14) is not well-behaved in the 
coordinates adapted to (2.15), i.e. the natural geodesic coordinates {p,yi), where yi are local 
coordinates on extended to coordinate functions on M to be invariant under the flow of Vp. 
Further, with this choice of conformal gauge, it is difficult to control the scalar curvature R of g. 

It is simpest and most natural to choose a conformal gauge of constant scalar curvature, (although 
other choices are possible). Thus, set 

(2.17) 9 = <7^9, 

where a is chosen to make R = const. In this gauge, the equation (2.13) for g simplifies to 

(2.18) n = {D*D)'^-'^Ric + L{D"g) = 0. 

The choice of constant for R is not important, but it simplifies matters if one chooses 

(2.19) R = R\x+ = = CO. 

77, — 1 
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The middle equality follows by taking the trace of (2.6), and combining this with the Raychaudhuri 
equation on g and (2.28) below. 

For the diffeomorphism gauge, we choose, as usual, harmonic coordinates Xa with respect to g; 

OXa = 0. 

It is assumed that the Cauchy data for xq are such that xq is a defining function for near , 
and we relabel xq = t so that the coordinates are {t,Xi), i = 1, ...,n. As usual, Greek letters are 
used for spacetimc indices, while Latin is used for spatial indices. Equivalently, but from a slightly 
different point of view, given arbitrary local coordinates x^, with xq a defining function for the 
boundary, the condition that Xa is harmonic with respect to g is 

(2.20) DXa = dar^ + ^r^Vdag^. = 0. 

In the coordinates Xq,, the metric g in (2.15) becomes 

(2.21) g = gooidt^) + 2goidtdxi + gijdxidxj, 
where 

(2.22) gQQ = -{dtpf, goi = dtpdip, gij = dipdjp + {gp)ij. 
Similarly, for g, one has 

(2.23) g = gooidt^) + 2goidtdxi + gijdxidxj, 
with gap = cr'^gap- 

As long as the coordinates are ^'-harmonic, the Ricci curvature has the form 

RicajB = -^g'"'dij,dugap + Qapig, dg), 

Similarly at leading order, the Laplacian D*D has the form g^^d^d^, in harmonic coordinates. Thus, 
with these choices of gauge for the conformal and diffeomorphism invar iance, the equation H = 
has the rather simple form 

(2.24) {r''d^d,)'^gap + L{D^g) = 0. 

This is an iV X system of PDE's for g^is which is diagonal, i.e. uncoupled, at leading order, 
N = (n+ l)(n + 2)/2. These choices for the conformal and diffeomorphism gauges are the simplest; 
however, they are not necessary and other choices, for instance gauges determined by fixed gauge 
source functions, cf. [12], could also be used. 

Having discussed the equations for the metric, we have left to determine the equations for p in 
(2.15) and a in (2.17). The fact that p is a geodesic defining function for g, i.e. \Vp\g = —1, implies 
that 

dtig'^^'dapOpp) = 0, 

or equivalently, 

(2.25) dt{a^r''dapdf3p) = 0. 

To derive the equation for a, the equation for the Ricci curvature relating g and g is 

(2.26) Ric = Ric+{n-l)^-^ + {—-n\dloga\'^}g. 

a a 

Taking the trace gives the equation relating the scalar curvatures as 

Da 
a 



(2.27) a-^R = R + 2n— - n(n + l)|£i log C7||. 
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Using the formula analogous to (2.26) relating the Ricci curvature of g and g, together with the 
fact that g satisfies (2.1) and p is a geodesic defining function gives 

(2.28) R = -2n^ = -2nRic{T, T), 

where T = dp = —Vgp. (Observe that the middle term in (2.28) is degenerate at T"*", since p = 
there; however, the last term in (2.28) is non-degenerate at p = 0). Substituting (2.28) in (2.27) 
and using (2.26) gives then the equation 

^ + n(n- 1) Id log o-|| = -2na-^[Ric(T,T) + (n - 1)^-^(T,T)], 

or equivalently, 

(2.29) TT{a) - {Va^VrT) = ~-^aRic{T,T) - ^^^^'_^^ co - ^a|dloga||, 

where we have also used (2.19). 

The equations (2.24), (2.25), and (2.29) represent a coupled system of evolution equations for 
the variables {gai3,P,(^) on a domain U in (R"'+^)+ with coordinates {t,Xi); the boundary doU = 
U r]{t = 0} corresponds to a portion of X+. Written out in more detail, these are: 

(2.30) {r''di,d,)'^9ap = Li{D^g)^p, 

(2.31) ^°dtdtp + 2{V'p){didtp) = L2{Dp, Da, Dg), 

(2.32) iy^fdtdtcj + 2{y%){V'p)didta + {V'p){V^ p)didja = L^{Da, D^g). 

Here D'^w denotes derivatives up to order k in the variable w and V"/? denotes the a-component of 

Vp, Vp = g°'^dfipda- The terms Lj are lower order terms. Observe that the system (2.30) for the 
metric ga/S is a closed sub-system, i.e. it does not involve p or a. Moreover, although the equations 
(2.31) and (2.32) for p and a are coupled to each other and to (2.30), the system (2.30)-(2.32) is 
uncoupled at leading order. 

Following common practice, we now reduce the system (2.30)- (2.32) to a system of 1®* order 
equations. There is not a unique way to do this, but we will discuss perhaps the simplest method, 
which uses pseudodifferential operators. As usual, the domain doU C is viewed as a domain in 
the n-torus T" and the variables (gap, <7, p) are extended to functions on 7 x T". 

Recall that a system of 1^* order evolution equations 

m 

(2.33) dtu = Bj{t, X, u)djU + c{t, x, u) 

i=i 

is symmetrizable in the sense of Lax, cf. [20], [21], if there is a smooth matrix valued function 
R{t, u, X, ^) on M X MP X T*{T^) \ 0, homogeneous of degree in ^, such that i? is a positive definite 
p X p matrix with R{t, u, x, ^) ^ Bj{t, u, x)^j self-adjoint, for each (t, u, x, ^). It is well-known [20], 
[21] that strictly hyperbolic systems of PDE, diagonal at leading order, are symmetrizable. A 
symmetrizer R is given by i? = X^-Pfe-P^' where Pk is the projection onto the k^^ eigenspace of the 
symbol ^ Bj [t, u, x)^j , 1 < k < p. 

Proposition 2.1. There is a reduction of the system (2.30)-(2.32) to a symmetrizable system of 
1^* order evolution equations on I x T"'. 
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Proof: Consider first the closed system (2.30) for g. This system is not strictly hyperbolic; the 
leading order symbol is diagonal and has two distinct real eigenvalues, each of multiplicity (n + l)/2. 

~ n+l 

However, the eigenspaces of the symbol of □ 2 vary smoothly and do not coalesce. Thus the 
operator □"is strongly hyperbolic, cf. [12] and references therein. In these circumstances, it is 

essentially standard that the operator □'^ is symmetrizable, for any k; for completeness we sketch 
the proof following [20, §5.3]. 

Let u = gap be the variable in M^, iV = (n + l){n + 2)/2. Write □, (we drop the tilde here and 
below), in the form 

1 

j=o 

where Aj is a differential operator in x, homogeneous of order 2 — j, depending smoothly u. Then 

n 

(2.34) [{g'')-'n](^+'y' = 0^' - Y.Bj{u,DM, 

j=0 

where Bj are differential operators in x, homogeneous of order n + 1 — j. Set Uj = 9^A"~-^n, for 
j = 0, ...,n, where A = (1 — A)^/^ and A is the standard Laplacian on T". Then (2.30) becomes 

n 

dtUj = Auj+i, 0<j<n, dtUn = ^ Bj{Pu, Dx)A^~'^Uj + C{Pu), 

j=0 

where Pu = D^u involves n up to n derivatives. More precisely, for /? + j < n, dxd{u = dxA^~"'Uj; 
for example A~"?xo = u. This is a system of l'''^ order pseudodifferential equations in the variables 
u = {uj}, j = 0, n of the form 

(2.35) dtu = L{Pu,Dx)u + C{Pu). 

The eigenvalues Xiy{w,^) of the matrix L{w,^) are the roots of the characteristic equation t'^+'^ — 
^ i?j(tf,^)T-^, (up to an overall factor of i). Hence, from (2.34), one sees that for each {w,$,), 
^ 7^ 0, there are two distinct roots, each of multiplicity (n + l)/2. The eigenvalues vary smoothly 
with (u),^) and remain a bounded distance apart on the sphere |^| = 1. The same is true of the 
corresponding eigenspaces. Hence, the system (2.35) has a symmetrizer R constructed in the same 
way as following (2.33), cf. [20, Prop.5.2.C] or [21, Prop.16.2.2]. 

Next we show that the equation (2.31) is also symmetrizable. Let (pi = —2\/^p/'(f^. Introducing 
the vector variable v = {p, po, pn) with po = dtp, Pj = djp, the equation is equivalent to the 
system 

dtp = Po, dtpo = ^ (pjdjpo + c{t, X, v), dtpj = djpo. 

This has the form 

n 

(2.36) dtv = ^'^Bj{x,t,v)djV + c{x,t,v), 

i=i 

where Bj is an (n + 2) x (n + 2) matrix with cpj in the (2,2) slot, 1 in the (j + 2,2) slot, and 
elsewhere. The system (2.36) is coupled at lower order to the equations (2.30) and (2.32) for g and 
a respectively, in that Bj depends on g to order 0, while c depends on g and cr to order 1; for the 
moment, these dependencies are placed in the (x, t) dependence of Bj and c. 

The matrix J2 ^jCj the entry J2 (Pj^j in the (2,2) slot, ^j in the (j + 2, 2) slot for 3 < j < n, 
and elsewhere. By a direct but uninteresting computation, it is straightforward to see that this 
matrix is symmetrizable in the sense following (2.33). 
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Essentially the same argument shows that the equation (2.32) for a is again symmctrizable. Thus 
let w = (cr, (To, ...,<7„) with ctq = dfU, aj = djCj. The equation (2.32) is equivalent to the system 

dta = ao, dtao = '^(f)jdjao + ^'(pijdjai + c{t,x,v), dtPj = djpo, 
where = -2(VV)/|V°p|, Vij = (VV)(V^p)/|V°|Op. This system has the form 

n 

(2.37) dfW = Bj(x, t, w)djW + c{x, t, w), 

i=i 

where Bj is the (n + 2) x (n + 2) matrix with cpj in the (2, 2) slot, 1 in the (j + 2, 2) slot, and tpij 
in the (i + 2,j + 2) slot. The system (2.37) is again coupled at lower order to the equations (2.30) 
and (2.31) for g and p respectively, in that c depends on g to order 2, while Bj depends on g to 
order and p to order 1. Again a straightforward but longer (uninteresting) computation shows 
that the matrix ^ Bj^j is symmetrizable in the sense of (2.33). 

One may then combine the three systems (2.35), (2.36), and (2.37) to a single large system 
in the variable U = {u,v,w). The resulting system is then a symmetrizable system of 1^* order 
pseudodifferential equations, cf. [20], [22]. ■ 

Next consider the Cauchy data for the system (2.30)-(2.32). If one is interested in general solu- 
tions of this system, then the Cauchy data are essentially arbitrary, subject only to the constraint 
equation H{Vp, •) = on X+. However, as will be seen in Proposition 2.3 below, it is the speci- 
fication of the Cauchy data which determines the class of conformally Einstein metrics among all 
solutions of (2.30)-(2.32). This is of course closely related to the FG expansion (2.8) of the dS 
metric g. 

The Cauchy data for a are 

(2.38) cr = 1, and dta = at X+, 
while the Cauchy data for p are 

(2.39) p = 0, and dtp = 1 at X+. 

For the metric g, the closed subsystem (2.30) is of order n -|- 1, so Cauchy data are specified by 
prescribing {dt)^gai3, k = 0,1, ...,n at Z~^. We compute the data inductively. First, the condition 
(2.38) implies that gij = gij at X+. Thus at order 0, set 

(2.40) goo = -1, 9oi = 0, gij = jij at X+, 
since pj = at X"*". 

At order, (2.38) and (2.39) together with (2.23) show that dtgij = dtgij at I+, and the FG 
expansion (2.8) gives dtgij = at X"*". Thus, set 

(2.41) dtgij = at X+. 

(This condition, and related ones below, are necessary to obtain Einstein metrics). The first 
derivatives of the mixed components goa of g are determined by the requirement that the coordinates 
Xa = {t, Xi) are harmonic at X+ with respect to Ij, i.e. for each j3, 

(2.42) d^r" + Ir^r^dagi.. = o. 

Via (2.40)- (2.41), this determines dtgoa at X+. 

At 2^^^^ order, the equation (2.27) implies, using the normalization (2.19), that dfa = at X"^. 
Also, dt{piPj) = 0, and hence, from the FG expansion (2.8), we set 

(2.43) d'^gij = 2ff(2) at X+. 
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The 2 derivatives dtgoa arc then determined by (2.43), the lower order Cauchy data, (2.38)-(2.41), 
and the t-derivative of (2.42) at t = 0. 

At S''^ order, suppose first n = 3, so that diniM = 4. Then a straightforward computation, using 
(2.8), the Raychaudhuri equation on g and (2.28), shows that dtR = Gntrg^^^ = at X+, where the 
last equality follows from (2.7). Hence, (2.27) gives dfa = 0. Similarly, (2.39) gives df{piPj) = 0. 
Thus, set 

(2.44) dfgij = 65(3) at X+. 

This term is free or unconstrained, subject to the transverse-traceless constraint (2.7). As before 
the mixed term at order 3, dfgoa is determined by taking two t-derivatives of (2.42) at t = 0, and 
using (2.44) together with the determination of the lower order Cauchy data. 

Suppose instead n > 3 and hence n > 5. Then g^^-^ = and same arguments as above give 

(2.45) dfgij = at T+, 

with dfgoa again determined from two t-derivatives of (2.42), (2.45) and lower order Cauchy data. 

At 4*^^ order on Z"*", (assuming n > 5), by (2.27) and the fact that d^a = 0, A; < 3 on X^, one 
has dfa = —2nd^R, and again computations as above then give d^R = 24ntrg(^4y Also, taking 
i-derivatives of (2.31) or (2.25) and using (2.38)-(2.39) shows that df{piPj) = at 1+ . It follows 
that, at t = 0, 

(2.46) dtgij = 24(5(4) - 2nHrg^^))^. 

Again, diTj^a is determined by taking three t-derivatives of (2.42) at f = 0, and using (2.46) with 
the determination of the lower order Cauchy data. 

At order, suppose n = 5. As in the case n = 3, dfR = ctrg{^-) = at X"*" while df{piPj) = 0. 
Hence, as in the case n = 3, 

(2.47) 5,% = (5!)<7(5), 

which is freely specifiable, subject to the transverse-traceless constraint. As before the mixed term 
at order 5, dfgoa is determined by taking four t-derivatives of (2.42) at t = 0, and using (2.47) with 
the determination of the lower order Cauchy data. 
If n > 5, then as before, 

(2.48) dfgij = 0, 

with Ofgoa again determined from (2.42). It is clear that one can continue inductively in this 
way to determine the Cauchy data d^ga/3 up to order n. Since df{piPj) 7^ at i = 0, these and 
higher derivative terms contribute to the Cauchy data at order 6 and above. However, one sees 
by differentiations of (2.31) and (2.42) that these terms are all determined by lower order Cauchy 
data for g. 

In sum, Cauchy data for gap, i.e. the data d^'gap, < k < n, are determined by the Cauchy data 
(2.38), (2.39) for p and cr, the equations (2.31)-(2.32), (or (2.25), (2.27)) for p and a, the harmonic 
equation (2.42), and the coefficients (7(fc) in the FG expansion (2.8). Thus, the Cauchy data are 
uniquely determined in terms of the free data 

(2.49) 7 = 5(0), and r = 5(„), 

which are arbitrary, subject to the transverse-traceless constraint (2.7) on g(^n) and the constant 
scalar curvature constraint (2.16) or (2.19) on the representative 7 G [7]. Abusing notation slightly, 
we will call (7,r) Cauchy data for g^p, since this data determines the rest of the Cauchy data 
Ojtdap, < k < n uniquely. 

The analysis above then gives: 
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Proposition 2.2. The system (2.30)-(2.32) for {gap,p,a) is well-posed in H^+'^{T+) x H%I+), 
s > f + 2. Thus, given Cauchy data (7,r) G {H'+'\I+), H'{I+)) satisfying (2.38), (2.39) and 

(2.49) , and satisfying the constraints (2.7), (2.19), there is a unique solution (gap-, Pif^) of (2.30)- 
(2.32) with 

(2.50) (5?a/3, P, a) G C(/, H'^^{I^)) n C-(/, H\I^)). 

Further, if^jg^Ts) is a continuous curve in H^~^'^{I^) x H^(I^), then the solutions (gs,Ps,o's) vary 
continuously with s. 

Proof: In the local coordinates {t, Xi), the system (2.30)-(2.32) is symmetrizable and so for given 
Cauchy data on T", it has a unique solution on / x satisfying (2.50), with T" in place of X"*", 
cf. [20, §5.2-5.3] or also [21], [22]. The existence of such local solutions holds for s > § + 1. By 
restriction, one thus obtains local solutions on the domain / x C / x T", for [/ C as preceding 
(2.30). 

To prove these local solutions obtained from domains on I x patch together to give a unique 
solution on I x X"*", it is necessary and sufficient to prove that the system (2.30)-(2.32) has finite 
domains of dependence (or equivalently uniqueness in the local Cauchy problem). This is well- 
known to be true for symmetric hyperbolic systems of PDE's, (cf. [22, §IV.4] for example); however, 
the reduction of (2.30) is to a symmetric system of first order pseudo differential equations, for which 
finite propagation speed is not true in general. Nevertheless, standard methods show that solutions 
of (2.30)-(2.32) do have local uniqueness in the Cauchy problem. 

Thus, consider the closed subsystem (2.30) for g. A standard argument using the mean-value 
theorem shows that it suffices to prove local uniqueness in the Cauchy problem for the associated 
linear equation L(g) = 0, where the coefficients of (2.30) are frozen at a given metric cf. [12] for 
instance. The linear operator L is self-adjoint at leading order, and hence by Proposition 2.1, one 
has existence and uniqueness of solutions to the Cauchy problem on / x T" for the adjoint equation 
L*u = (p. It is then well-known, cf. [22, Thm.IV.4.3], that this suffices to prove local uniqueness 
in the Cauchy problem for L, and hence for the system (2.30). Given the local uniqueness for g, 
exactly the same method can be applied to (2.31)-(2.32) to give local uniqueness for p and a. 

It follows that the system (2.30)-(2.32) does have finite propagation speed, and hence the local 
solutions patch together uniquely to give a unique global solution on 7 x satisfying (2.50). As is 
well-known in the case of Einstein metrics, (cf. [12] for example), the patching of local coordinate 
charts requires an extra derivative; the same analysis holds here, so that we assume s > ^ -\- 2. 

A lower bound for the time of existence / = [0, io) of the solution depends only on an upper 
bound on the norm of the initial data in x H^. This implies that the last statement is an 

immediate consequence of the existence and uniqueness theorem. 

■ 

Given a solution (ga/s, p, c) of the Cauchy problem, one may then construct the "physical" metric 
g by setting 

(2.51) g = p-^a-^g. 

Since a is bounded away from and 00 near X"*" and p is a geodesic defining function, it is easy 
to see that the metric g is future geodesically complete, i.e. geodesically complete to the future 
of some Cauchy surface E C (M"+^,5). However, it is not so clear that the metric g is Einstein, 
or equivalently that the metric g is conformally Einstein. For this, one needs to verify first that 
the gauge condition (2.42) that the coordinates remain harmonic, is preserved for the solution g. 
If this is so, it then follows that ^ is a solution of the equation (2.11). Secondly, the equation 
(2.11) admits solutions which are not conformally Einstein, and so one needs to verify that the 
constructed solution g actually is conformally Einstein. 
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Both of these conditions on g can be verified by computation; however, the computations will 
be somewhat long and involved. Instead, we verify both conditions together by using the following 
simple conceptual technique, based on analyticity. 

Proposition 2.3. Any solution (gafi-, Pi^^) of the Cauchy problem in Proposition 2.2 defines an 
Einstein metric in dS'^ via (2.51). 

Proof: Suppose first that the free data (7, ^(n)) are analytic on 2r+ (or on a domain in X+). 
As noted above, the expansion (2.5) then converges to Qp and gives a solution, denoted qe, to the 
Einstein equations (1.1); this metric has the form (2.4), with compactification qe., in a neighborhood 
of . In particular, both qe and p are analytic near . Moreover, since the coefficients of the 
equation (2.27) defining a are then also analytic, and since the Cauchy data (2.38) for a are analytic, 
the Cauchy-Kowalewsky theorem ([21, §16.4]) shows that a and hence IjE = are analytic near 

. Of course the metric qe is a solution of (2.11), and hence a solution of (2.30)-(2.32). with 
Cauchy data determined by (2.38)-(2.39) and (7,(7(n)). 

On the other hand, let {jj, p, a) be the solution to the Cauchy problem (2.30)-(2.32) with Cauchy 
data determined by (2.38)-(2.39) and (7, </(„)) given by Proposition 2.2. Since (gE,P,cr) is also a 
solution of this Cauchy problem, with the same initial data, it follows from the uniqueness part 
of Proposition 2.2 that g^; = ^ in a neighborhood of I'^. Moreover, the results in [1] show that 
qe remains analytic within its globally hyperbolic development, so that gE = 9 everywhere in the 
domain of g. Hence g = p~^a~^g is an Einstein metric in dS~^, realizing the given data on (In 
particular, this shows that the harmonic gauge (2.42) must be preserved for analytic initial data). 

Now analytic data (7, 5(„)) are dense in x data with respect to the x topology 

on The Cauchy stability given by Proposition 2.2 implies that if {'Ji,g(n),i) is a sequence of 
analytic Cauchy data converging in H'^^^ x if* to if*"*"" x data (7, g(n)), then the corresponding 
solutions (gi,Pi, ai) also converge to the (unique) solution (g, p, a) with Cauchy data (7, (?(„)). Hence 
the metric g in (2.51) is Einstein. 

■ 

Propositions 2.2-2.3 give the existence of an Einstein metric {M^~^^,g) € dS~^ , with arbitrarily 
prescribed asymptotic behavior (7, r) on I+, subject to the constraints (2.7) and (2.19). Suppose 
[7'] = [7] and [r'] = [r] as following (1.5), so that 7' = A^7 and r' = /(A)r. Given A, there exists 
a diffeomorphism cp : M ^ M, with 4'\j+ = id, such that linip^Q (p* (p) / p = X~^, where p is the 
geodesic defining function determined by 7. Setting g' = (j)*g, one has g' = p^(l)*{g) = 'P*{9){ ^f{p) )^- 
Hence, the boundary metric of g' equals 7'. Similarly, by the uniqueness, any Einstein metric 
(M^+^,5') with Cauchy data (7',r') satisfying (1.5) differs from (M"+^,5f) by a diffeomorphism cf) 
equal to the identity on 2"^. This completes the proof of Theorem 1.1. 

■ 

Remark 2.4. Theorem 1.1 is formulated as a global result, in the sense that future conformal 
infinity X+ is a compact smooth manifold. However, Propositions 2.1-2.3 all hold locally, by the 
finite propagation speed of the system (2.30)-(2.32). Hence, Theorem 1.1 also holds locally, where 
X+ is an open manifold with a finite number of local charts. Of course the uniqueness statement 
then holds only within the domain of dependence of the initial data. 

Many of the standard solutions of the Einstein equations (2.1) have Z non-compact; this is the 
case for instance for the dS Schwarzschild metrics. 

Proof of Theorem 1.2. 

Let (M, 50) be a dS Einstein metric in dS"^ with Cauchy data ([7"^], b^^)]) and ([7~], [5(^)1) 
induced on X"*" and T" respectively. Thus, there exists a smooth defining function : M ^ R such 
that ^0 = ^^90 extends to a metric on M = M U X+ U X~ ~ / x S; here E is a Cauchy surface for 
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(M, 5'o) and / is a compact time interval. We will choose to be a geodesic defining function p in 
a neighborhood of X+ and X" so that go is C{I x n C"(/ x up to M. The choice 

of n defines representatives (7+,5f+p, (7",5(~)) in the conformal classes ([7+], b^)]), ([7"], [%)])• 
In the following, we work in the x topology. 

Let be an open neighborhood of Cauchy data on X+ containing the given data {l'^,9(n))- 
Then for all data {'y~^,g^^^) G W"*", there exists T < 00, depending on U'^, such that the maximal 
globally hyperbolic dS Einstein metric {M"'~^^,g) having Cauchy data {'1~^t9^^) given by Theorem 
1.1 is defined on [T, 00) x E; here the time factor is proper time t = — log ^, where p is the geodesic 
defining function. The Cauchy data of such solutions ^ at S = {T} x S then forms an open set 
in the space of Cauchy data for the Einstein equations (2.1) on S. By passing to an open 
subset C if necessary, the Cauchy stability theorem for the (standard) Einstein equations 
implies that the maximal globally hyperbolic development of any g with data in contains the 
region [— T, T] x E, and induces again an open set of Cauchy data at {— T} x E. Then, as 
above with there is an open set of Cauchy data on X^ whose future development gives a 
non-empty open subset of V^^. Combining these three unique developments gives a global solution 
(M"+^,y) € dS^, which completes the proof of Theorem 1.2. 

■ 

Remark 2.5. The proof above shows that the space dS"^ is also stable with respect to perturbations 
of the Cauchy data (E,7, i^T), (satisfying the constraint equations), on a compact Cauchy surface 
E C (M"+^9) in the x H^+^-^ topology on E. It is well-known, cf. [4] that this is not the 

case for perturbations of asymptotically flat Cauchy data when A = 0, in that smoothness of the 
resulting space-time at conformal infinity is lost for generic perturbations. 

Proof of Theorem 1.3. 

This result is proved forn = 3 in [2], and it is pointed out there that the same proof holds provided 
one has the Cauchy stability result of Theorem 1.1, (i.e. Priedrich's result [8] when n = 3). Given 
then Theorem 1.1, the proof of Theorem 1.3 is exactly the same as that given in [2], to which we 
refer for details. 



Remark 2.6. This paper has focussed on the de Sitter case A > mainly for simplicity, but also 
because there are no direct analogues of Theorems 1.2 or 1.3 when A = or A < 0, due to the more 
complicated nature of conformal infinity. Nevertheless, one expects that the analogues of Theorem 
1.1 for A = and A < 0, as formulated and proved by Friedrich [8] in the case n = 3, hold for 
all even dimensions. When A < 0, future space-like infinity X"*" is replaced by time-like infinity X, 
while when A = 0, is replaced by future null infinity. 
We hope to discuss these situations elsewhere. 

Remark 2.7. We close with a brief remark on the applicability of the methods used above to the 
Einstein equations coupled to other matter fields. As noted above, up to multiplicative constants, 
the tensor H is the metric variation, or stress-energy tensor, of the conformal anomaly or of the 
Q-curvature. The conformal invariance of these functionals corresponds to the conformal invariance 
of H. For functionals containing the metric coupled to other fields which are conformally invariant, 
and whose field equations are symmetric hyperbolic, it seems very likely that the methods used 
above will again lead to a well-posed Cauchy problem at X"*", as in Theorem 1.1. In dimensions 
3-1-1, this is the case for the Einstein equations coupled to gauge fields, i.e. Einstein-Maxwell or 
Einstein- Yang-Mills fields. Theorem 1.1 has already been proved in this situation by Friedrich [11], 
and so at best one would have a different method of proof of this result. In higher dimensions, the 
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EM or YM action is not conformally invariant, and it is less clear if the method can be adapted to 
this situation. 
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